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Abstract 


Laminate and structural mechanics for the analysis of laminated composite plate structures with 
piezoelectric actuators and sensors are presented. The theories implement layerwise representations 
of displacements and electric potential, and can model both the global and local electromechanical 
response of smart composite laminates. Finite-element formulations are developed for the quasi-static 
and dynamic analysis of smart composite structures containing piezoelectric layers. Comparisons with 
an exact solution illustrate the accuracy, robustness and capability of the developed mechanics to 
capture the global and local response of thin and/or thick laminated piezoelectric plates. Additional 
correlations and numerical applications demonstrate the unique capabilities of the mechanics in 
analyzing the static and free-vibration response of composite plates with distributed piezoelectric 
actuators and sensors. 
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1. Introduction 


The development of a new class of "smart" composite materials and adaptive structures with 
sensory/active capabilities may further improve the performance and reliability of aeronautical 
structural systems. Such materials will combine the superior mechanical properties of composite 
materials, as well as, incorporate the additional inherent capability to sense and adapt their static and 
dynamic response (adaptive structures), or continuously monitor the type, location, and extent of 
eminent damage (health monitoring). However, this effort requires the development of admissible 
mechanics entailing capabilities to model the unified electromechanical response of sensory/active 
structures including the coupling between sensors and actuators. Additionally, the mechanics should 
address local through-the-thickness effects, such as the evolution of complicated stress-strain fields 
in smart composites and interfacial phenomena between the embedded micro-devices and passive 
composite plies in a smart laminate. The present paper presents the development of layerwise 
mechanics for the dynamic analysis of smart composite plate structures with embedded piezoelectric 
sensors and actuators. 

There have been many theories and models proposed for the analysis of laminated composite plates 
co ntainin g active and passive piezoelectric layers (see e.g. Lee and Moon (1989), Lee (1990), Lee 
and Moon (1990), Wang and Rogers (1991), Crawley and Lazarus (1989), Tzou and Tseng (1990), 
Lammering (1991)). Most of these theories use simplifying approximations attempting to replicate 
the induced strain or electric fields generated by a piezoelectric layer under an external electric field 
or applied load. An exact solution by Heyliger and Saravanos (1995) for piezoelectric laminated 
plates has shown that the electric and elastic field distributions are often poorly modeled using 
amplified theories. Variational methods and finite element models for piezoelectric solids have also 
been reported by Allik and Hughes (1970), Naillon et al. (1983) and Ha et al. (1992). Related work 
has been reported for infinite piezoelectric laminates by Pauley (1974), and for finite elastic laminated 
beams and plates with induced strain actuation by Robbins and Reddy (1991, 1993). The work in this 
paper builds on previous layerwise laminate theories (Reddy (1987), Robbins and Reddy (1991), 
Robbins and Reddy (1993)), as well as, on that of the authors (Saravanos (1993), Saravanos (1994), 
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Heyliger et al. (1994), Saravanos and Heyliger (1995)) to consider the complete dynamic electr- 
omechanical response of smart piezoelectric plate structures under external mechanical or electrical 
loading. 

The paper describes layerwise laminate theories for composite structures containing embedded 
piezoelectric layers as sensors or actuators. According to these laminate theories, which will be called 
interchangeably in the remaining paper as layerwise or discrete-layer theories, each layer is modeled 
using independent approximations for the in-plane displacement components and the electrostatic 
potential in a unified representation, as mandated by the linear theory of piezoelectricity. The laminate 
theories assume either constant or variable transverse displacements through the laminate thickness. 
Finite element solutions for the dynamic response of smart structures are further developed. 
Application studies verify the accuracy and demonstrate the advantages and versatility of the 
mechanics in analyzing the dynamic response of composite plates with sensory and active 
piezoelectric elements. 


2. Piezoelectric Laminates 

This section describes the developed mechanics for composite laminates with embedded piezoelectric 
sensors and actuators. The coupled material equations for each ply are first presented in a unified way 
which may represent either piezoelectric or passive composite layers. Discrete-layer approximations 
are subsequently defined and the resultant piezoelectric laminate theories are presented. A standard 
laminate coordinate system is assumed, such that, x and y are the in-plane axes and the z-direction 
coincides with the thickness dimension of the laminate. The general problem considered in this paper 
is to determine the dynamic behavior of the elastic and electric field components throughout the 
laminate under an applied mechanical or electrical loading. Each ply of the laminate can be composed 
of a purely elastic, piezoelectric, or conducting material. The forcing function may be a combination 
of applied tractions, electric potential, surface displacements, or charge. 
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Material Representation 


The constitutive equations of a piezoelectric material are given by Tiersten (1969), 

S , ■ * d a£k 

D i - fyj * 


or equivalently, 


°i ' Cy Sj - e aFk 
D t - ejSj ♦ t 0 E k 


( 2 ) 


where: ij = and k,l= 1,...,3; o t and S ; are the mechanical stresses and engineering strains in 
vectorial notation; E* is the electric field vector; D, is the electric displacement vector; s s and C a are 
the elastic compliance and stiffness tensors; d a and e^ are different forms of the piezoelectric tensor, 
and e* is the electric permittivity tensor of the material. Superscripts E, o, and S indicate constant 
electric field, stress and strain conditions, respectively. The electric field vector Ej. is the gradient of 
the electric potential 4> 

E k - -di>/dx k (3) 

The materials are assumed to be monoclinic class 2 crystals with a diad axis parallel to the z axis. The 
poling direction of the material is assumed coincident with the z axis. The assumed material class is 
general enough, such that eqs. (1) or (2) may encompass the behavior of off-axis homogenized 
fibrous piezoelectric plies, as well as, passive composite plies [e]=0. 


The energy stored (electric enthalpy) in the piezoelectric layer includes the components of elastic 
strain energy, piezoelectric energy and electric energy (Tiersten (1969)), 

H t - ^.C^j - 


( 4 ) 



Discrete-Layer Laminate Theory 


Two new layerwise laminate theories (also called discrete-layer theories) are proposed in this paper 
using piecewise continuous approximations along the z-axis for both displacement and electric 
potential fields. Previous work has demonstrated the advantages of layerwise approaches in capturing 
local intrala minar and interlaminar effects in elastic (Reddy (1987)) and elastodynamic problems 
(Saravanos (1993), Saravanos (1994)) of composite plates. Consequently, the suitability of discrete- 
layer theories to represent the additional heterogeneity induced by the presence of embedded 
piezoelectric sensors/actuators in composite laminates can not be understated. 

In the present approach, the electric potential is included in the state variables. The laminate is 
subdivided in N-l sublaminates (or discrete layers), and continuous fields are assumed in each 
sublaminate, such that, piecewise continuous representations result through the thickness of the 
laminate (see Fig. 1) of the following form: 

N 

<xy^y £ u'(xyJ)W J (z) 

>1 

N 

v(xyjjy 52 v-fc yJ)V J (z) 

£ ( 5 ) 

w(xyjjy " J (xy,t)'P **(z) 

>1 

Jf 

4>(xyj,ty 52 ^{xyjtyf\z) 

>1 

where u, v, w are the displacements along the x, y and z axes; superscript) indicates the points z* at 
the beginning and end of each discrete layer, u*, v*, w*, are the respective displacements and electric 
potential at each point z* (see Fig. 1); Y j (z) are interpolation functions. Two unique advantages of 
the method are obvious: (1) the complete electromechanical state of the smart laminate is represented; 
and (2) the formulation entails the inherent option to select the detail of approximation in both electric 
and displacement fields. At the lower limit of one discrete-layer (N=2) the method reduces to "single- 
layer" type of laminate theories, and for linear Y(z) it reduces to the first-order shear theory. 
Moreover, the through the thickness displacement can be assumed constant (Y w (z)=l) or variable, 
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which establis he s the difference between the two laminate theories. The first, which will be identified 
as the w-constant theoiy, results by assuming Y w (z)=l in eqs. (5) and is effectively a plate laminate 
theory. The second, which is identified as the w-variable theory, is a three-dimensional laminate 
theory. In the present paper only linear interpolation functions Y(z) were considered. 

In the context of eqs. (5), the engineering strain and electric field components become, 

S{xyj,ty E s/(xyj)V J (z) .- 1 , 2,6 

>1 

s^x^y E £ 3 Wo*7(*) 

>1 

s/ixy^ty E slixy^/zy^xyjye^z) .- 4 ,s (6) 

>i 

N 

E/ix^ty E E i i (xyj)'i J {z) i-ia 

>1 
N 

E 3 (xyj,ty E E k x ’y* t )' f 
>1 

where the {S’} and (E’} are respectively the generalized strain and electric field vectors defined as 
follows: 

sf-u* sj-v' si-u^vj, si-w 1 sl-v * sj-u J 

s{-w' (7) 

e{.-^ E*-** 


The electric enthalpy of the laminate is by definition, 

£ *i * (8) 

where h is the laminate thickness. By combining eqs. (4,6,8) and integrating through-the-thickness, 
the energy stored in a piezoelectric laminate is obtained as a quadratic expression of the generalized 
strain/electric field and the generalized laminate stiffness, piezoelectric and dielectric permittivity 
matrices, 
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( 9 ) 


N N 
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In the above equation. A"”, B"“, and D™ are the generalized laminate stiffiiess matrices, 
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E™ are the piezoelectric laminate matrices. 


iV-4.5 


ET. E f‘ u e'f m (z)W’ , /z)dz «-3 > 1 , 2,6 
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(ii) 


and G"“ are the laminate matrices of electric permittivity, 

L 

Gr. E f ^e tf T"(*)¥"(x)dr 1-7-U 

M 

Z, 
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where L is the number of plies in the laminate. 


( 12 ) 
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3. Piezoelectric Structures 


The starting point for the equations of motion of a composite structure with embedded active/sensory 
piezoelectric layers in a variational form, is the Hamilton's principle for a piezoelectric medium 
(Tiersten (1969)) 


if ^pu/i r H(.S u £J 6 *>S-0 


(13) 


Here t is time, V and S are the volume and surface occupied by and bounding the solid, overbar o 
and overbar q are the surface tractions and surface charge, respectively, 6 is the variational operator, 
the overdot represents differentiation with respect to time, and H, represents the local electric 
enthalpy. 

Integration through-the-thickness yields the Hamilton's principle (equations of motion) in terms of 
the generalized laminate quantities defined in the previous section 

N }? 

6 f ‘dtf (-£ £ ‘dtf (5 k t>u i -qt,*)dS-0 (14) 

» *A 2 M >1 Jft J s 


where, A is the mid-plane area occupied by the plate, H L is the electric enthalpy of the laminate given 
by eq. (9), and P* j are the generalized densities (per unit area) of the laminate given by 

P*- £ f* (15) 


To develop structural solutions, local in-plane approximations of the generalized electromechanical 
state (displacements and electric potential) in eq. (5) are proposed of the following type. 
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( 16 ) 


where 17 is the value of the generalized displacement component LP corresponding to the i-th in-plane 
interpolation function Rj(x,y), and so forth. This formulation forms the basis of the finite element 
technique, which is subsequently described. 


Finite Element Formulation 


For structural problems with general boundary, geometry and material configurations, the local in- 
plane approximations in eq. (16) are used to develop finite dement based solutions. Substituting these 
approximations into the equation of motion (14) and collecting the coefficients, allows the governing 
dynamic equations of the structure to be expressed in a discrete matrix form as, 

[A/ u ] [0] [0] [0] 

[0] \M n ] [0] [0] 

[0] [0] [WJ [0] 

. [0] [0] [0] [0] 




[r„] [/r 12 ] [* 13 ] [r 14 ] 

{U) 


W) 

tn 

► * 

[* 2 ,J [r j [**] [rg 

IV) 



{W} 

♦ 

[* 31 ] [* 3 J [KJ 

{W} 


(7,(0) 



[*<,] t r «l 

{*> 


.((2(0} 


The dements of these matrices are calculated from the generalized laminate matrices defined in eqs. 
(10-12) as determined by the variational statement. The nature of the submatrices depends also on 
the approximation used for w(z). For variable w, the structure of [K* 3 ] is similar to those of the other 
matrices. For constant w, the submatrices within [K 13 ], [K 23 ], and [K 43 ] are column vectors and those 
in [K 33 ] become scalars. In general, the submatrices within each [K ij ] are each of order (N), while the 
matrices [K* j ] depend on the order of the in-plane approximation (M). 
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Based on the above formulation, a 4-node finite element was developed with linear in-plane shape 
functions. Reduced integration (single point Gauss quadrature) was implemented in the calculation 
of the shear stiflness terms to eliminate overstiffening at low thickness. The final representation of 
the coupled dynamic system can be expressed in fairly compact form: 

[WJ o J {i?> I 1*J I J {U) | | im) | 

Assuming that both sensory and active piezoelectric layers are embedded into the structure, the 
electric potential vector is subdivided in a free or sensory component <t> F representing the voltage 
output at the sensors, and a forced or active component 4> A representing the voltage imposed on the 
active piezoelectric layers, such that {<$*} = { 4> F j c|> A }. Separating the active and sensory potential 
components in eq. (18), the equations take the following form 

p/J olj {«?>! firj p*f {«> lj (19) 

. o D.IW [*£i ix%) 1 f+ F >J {Q f V))-[k%\{* a } 

where superscripts F and A indicate the partitioned submatrices in accordance with the selected 
sensory and active configuration, respectively. The left-hand side includes the unknown 
electromechanical response of the structure {u, <|> F }, that is, the resultant displacements and voltage 
at the sensors. The right-hand includes the excitation of the structure in terms of mechanical loads 
and applied voltages on the actuators. The electric charge at the sensors Q F (t) remains constant with 
time (practically open-circuit conditions) and is assumed equal to zero. 

Among the obvious advantages is the capability of the mechanics to model the response of the 
piezoelectric structure either: in "active" mode, that is, with specified voltages A<J> A applied across 
the piezoelectric layers to induce a desirable deflection/strain state; or in "sensory" mode where 
displacements or mechanical loads are applied to the structure and the resultant voltage or charge is 
monitored; or in combined "active/sensory" mode. 

Additional manipulation of eqs. (19) results in the following uncoupled dynamic equations for the 
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structural displacements and sensory voltages respectively 

[wjw ♦ a*j-[0**r l [0><*> • ^(owiOiO-'iO-tO^) (20) 

{4> f } - (21) 

Further inspection of eq. (20) reveals that the "induced-strain" approaches, in the presence of sensors, 
neglect the coupling effects on both the stiffness and the induced piezoelectric force. The above 
dynamic system may be solved to obtain, either the modal characteristics (free-vibration), or the 
forced frequency response or the transient response of the piezoelectric plate. 


4. Evaluations and Discussion 

Application cases on laminated composite plates with active and sensory piezoceramic layers are 
presented in this section. The free-vibration response of a simply supported cross-ply plate with 
surface bonded continuous piezoceramic layers is analyzed first. Additional evaluations are presented 
on the active quasi-static response and free-vibration response of laminated cantilever plates with 
distributed discrete piezoelectric patches. When applicable, correlations with published data are 
presented. In all cases, the standard laminate notation is augmented to indicate the lamination and the 
location of the piezoelectric material through the thickness, with the letter p indicating the 
piezoelectric layer. 

Hybrid Simply-Supported Plates. 

The free-vibration response of a 5-ply [p/0/90/0/p] laminated simply-supported square plate 
incorporating piezoelectric layers and composite plies was analyzed. Two thickness aspect ratios are 
considered corresponding to a thin square plate (a/h= 50) and a thick plate (<x/h=4), where a is the 
free span and h is the plate thickness, to illustrate the range of application of the mechanics. Results 
obtained with the finite-element models are compared with an exact piezoelasticity solution reported 
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by Heyliger and Saravanos (1995). The laminate configuration consists of a [0/90/0] Gr/Epoxy 
cross-ply sublaminate with composite plies each 0.267h thick, where h is the plate thickness. Two 
continuous PZT-4 layers of thickness 0. lh each are also bonded to the upper and lower surfaces of 
the laminate. The elastic, piezoelectric and dielectric properties of the laminate materials are shown 
in Table 1 . To comply with the reported results of the exact solution, all layers were assumed to have 
equal density (p=l kg/m 3 ). The imposed simple support conditions were: w i (x=0, y)= w i (x=a, 
y)=w»(x, y=0)=w i (x, y=a)=0, v i (x=0, y)= v>(x= a, y)=0, u^x, y=0)=u>(x, y=a)=0. In addition, the 
outer surfaces of the piezoelectric layers were forced to remain always grounded. Based on this, two 
sets of electric boundary conditions were considered for the inner surfaces of the piezoceramic layers: 
a closed-circuit condition termed with (C), with the potential forced to remain zero (grounded); and 
an open-circuit condition termed with (O), where the electric potential remains free (zero electric 
displacements). 

Fundamental Natural Frequency. The predicted fundamental frequency fj of the plate with various 
uniform mesh densities and various numbers of discrete layers through the thickness of the laminate 
is shown in Table 2. In the case of three discrete-layers (N=4), one discrete layer was used for each 
piezoelectricity and one for the whole composite sublaminate; in the case of 7 discrete layers (N=8), 
2 discrete layers were used for each piezoelectric ply and 1 for each composite ply; and in the case 
of 20 discrete layers (N=21), 4 layers were used for each piezoelectric and composite ply. For both 
thickness ratios, the predicted natural frequencies consistently converge above and below the values 
of the exact solution depending on the type of electric boundary conditions, which shows that the 
element can model thick piezoelectric composite structures, while it does not lock at low thicknesses. 
Although the agreement is very good, the predicted range of the piezoelectric effect on the 
fundamental frequency (the difference between (O) and (C) conditions) is much higher in the FE 
results than in the exact solution, and the reason for this difference remains a subject of investigation. 
The consideration of minimal discrete-layers through the thickness seems to provide reasonable 
accuracy in the calculation of frequencies for both thickness ratios. Yet, higher number of discrete- 
layers may be required in thick laminates to capture the local laminate response, as shown in 
subsequent examples. 
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Electromechanical State. The through-the-thickness distributions of the modal displacements and 
the normalized electric potential of the fundamental mode of the plate at open-circuit conditions (O) 
are shown in Figs. 2-3 for a thin and thick plate, respectively. The transverse displacement w and 
electric potential correspond to the center of the plate (x=a/ 2 , y= a/ 2 ), while the in-plane 
displacements u and v are at the midspan of the edges, e. g. (x=a, y=a/ 2 ) and (x=a/ 2 , y= a), 
respectively. An 8 by 8 uniform mesh and 7 discrete layers through the thickness were used. In the 
case of the thin plate (a/h= 50, Fig. 2) both w-constant and w-variable laminate theories and 
corresponding finite elements (FE) yielded identical results with the exact solution. In the case of 
the thick plate (a/h= 4, Fig. 3), both laminate theories (w-constant and w-variable) and corresponding 
finite elements yielded excellent agreement with the exact solution. The w-constant theory slightly 
underpredicted the in-plane displacements, and did not capture the slight variation in w (mostly due 
to the dja piezoelectric coefficient) in the piezoelectric layers. The variation of the electric potential 
in the composite plies was caused by the in-plane electric permittivity £ u and £22 of the composite 
plies. Clearly, the assumptions of simplified theories (simplified linear displacement and electric fields, 
single layer assumptions) are only valid for low thicknesses (Fig. 2 ), while significant deviations and 
errors may be observed in thick plates (Fig. 3). 

Modal Stresses. The predicted modal stresses through the thickness of the laminate are shown in 
Figs 4-5, for the thin and thick plate, respectively. The stresses were calculated at the points of 
respective maximum values, that is, at the center (<j„ a y ), the comer (<y v ), and the mid-edges of the 
plate (Cy^ < 7 , 0 ). An 8 x 8 uniform finite element mesh was used with 20 discrete layers through the 
thickness. Calculated stresses with the w-constant finite element model are shown in Figs. 4-5, 
corresponding to the integration points and elements closer to the respective stress location. Again 
for both thickness ratios, the agreement of the finite element predictions to the exact solution is 
excellent. Fig. 5 shows clearly that in thick laminates, the stresses exhibit substantial nonlinearity and 
slope discontinuities between each ply, which can not be captured by simplified kinematic 
assump tions The previous results and comparisons have demonstrated the excellent accuracy and 
robustness of the mechanics and the finite elements, and their capability to capture both the global 
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and local through-the-thickness response of thick and thin piezoelectric plates. 

Effect of Electric Conditions. Figs. 6a-d illustrate the effect in the modal in-plane displacements, 
electric potential, and stresses of the thick plate, when different conditions are imposed on the electric 
potential at the inner surfaces of the piezoelectric layers. Specifically, the cases of forced (equal to 
zero) and free electric potential are shown in Fig. 6, identified as closed-circuit (C) and open-circuit 
(O) conditions respectively. As seen in Fig. 6a the voltage conditions have a definite effect on the 
modal in-plane displacements, as a result of the different electric fields in the piezoelectric layers (see 
Fig. 6b). It is interesting to note that electric fields exist in the piezoelectric layers even with the 
closed-circuit conditions. The effects of electric conditions are more profound in the in-plane stress 
ct x (Fig. 6c) and the interlaminar shear stresses (Fig. 6d), and are attributed primarily to the 
differences in the piezoelectric stress component (see eq. (2)). 

Composite Plates with Discrete Piezoelectric Patches 

The present applications illustrate the natural capability of the finite element to model more realistic 
configurations of smart structures, that is: (a) discrete piezoelectric actuators/sensors, and (b) 
continuous electrodes on the free surfaces which force the sensory electric potential to remain 
uniform. The plates considered in this example and the corresponding meshes are shown in Fig. 7. 
The first is a [0/±45], Gr/Epoxy cantilever plate with 15 piezoceramic patches at each side built and 
tested by Crawley and Lazarus (1989); the second is a candidate [P2/0/90/45/-45], T300/934 
Gr/Epoxy cantilever plate with 8 square PZT-4 patches at each ride with properties shown in Table 
1. Both plates were modeled using the w-constant theory with 3 discrete layers. 


[0/±45], Cantilever Plate. Fig. 8 shows the predicted transverse deflection of the 0.83 mm thick 
[0/±45], Gr/Epoxy (AS4/3501) plate shown in Fig. 7a with 15 piezoceramic patches bonded on each 
side, each 0.25 mm thick. The deflections are induced by an applied uniform electric field of 394 
V/mm, of opposite polarity at the upper and lower piezoelectric patches. Measured data were 
reported by Crawley and Lazar us (1989), together with some material properties (additional material 
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properties were backcalculated as required from related data in the same publication). The following 
normalized form was used for the deflections consistent with the way the experimental results were 
published, 

r, - vrjc, T 2 - T 3 - (wyW^e 

where, w* w b w 3 are the transverse displacement along the midline and the two edges respectively, 
and c is the width of the plate. The three normalized displacements respectively represent or 
approximate the axial bending deflection, the transverse bending curvature, and the twisting angle due 
to bending-twisting coupling. There is very good agreement between predicted and measured 
transverse deflection T 3 . The agreement between predicted and measured transverse curvature T 2 
and twisting angle T 3 is fair to good. The present predictions and correlation trends agree remarkably 
well with similar results reported by Ha et al. (1992), hence, the differences were attributed to 
uncertainties in the material properties and possible specimen imperfections. 

[PJ/0/90/45/-45], Cantilever Plate. Figs. 9-10 show the predicted first two electromechanical modes 
of this plate (Fig. 7b). All piezoceramic layers were configured to act as sensors, with continuous 
electrodes on each surface, which were represented with equality constraints imposed on the 
respective electric potential. There is a characteristic modal voltage distribution corresponding to the 
characteristic distortion of each mode. The first mode (Fig. 9) is almost a pure first bending mode in 
the axial direction. The second mode (Fig. 10) is the first torsional mode in the axial direction. Both 
figures indicate the non-apparent relationship between deflected shape and the corresponding sensory 
voltage pattern. The slight asymmetries in the sensory voltage patterns are caused by the flexure- 
twisting coupling of the laminate. Figs. 9-10 clearly illustrate the various static and/or dynamic 
deformations which may be actively induced, for example bending or twisting, by applying the proper 
pattern of voltage on the piezoelectric elements. Note also that undesirable deformations, caused 
by laminate coupling, may be annihilated in this manner. The natural calculation of the characteristic 
voltage distributions is an obvious advantage of the method, because it provides the voltage patterns 
required to induce the corresponding deflection, or to construct the modal filter which can be used 
to isolate, sense and ultimately control the respective mode. 
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Overall, all previous application cases have demonstrated the accuracy, quality and versatility of the 
developed laminate theory and corresponding finite elements, and their capability to model a broad 
array of global and local responses in most practical configurations of active/sensory piezoelectric 
composite plates. 

5. Summary 

Laminate and structural mechanics and the corresponding coupled electromechanical models for the 
dynamic analysis of smart composite plate structures with embedded piezoelectric actuators and 
sensors were developed and described. The described laminate theories entail mixed displacement and 
electric field formulations and implement layerwise representations of the displacements and electric 
potential. A robust linear plate finite element, in terms of thickness, was also formulated. The 
mechanics have capabilities to simulate both sensory and active dynamic response of smart composite 
structures either at the global structural or the local laminate level. 

Correlations of predicted results with exact solutions of simply-supported piezoelectric composite 
plates have demonstrated the accuracy, robustness and versatility of the developed mechanics. The 
comparisons verified the capability of the present mechanics to model the free-vibration of thick 
and/or thin piezoelectric composite plates. Additional application studies have demonstrated the 
capability of the mechanics to model the active and sensory, static and dynamic response of smart 
composite plates with continuous or discrete piezoceramic devices. Finally, the applications have 
demonstrated the possibilities for dynamic displacement management provided by smart piezoelectric 
composite structures. 
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Table 1 Mechanical Properties 
(e 0 =8.85 1 O ' 12 farad/m, electric permittivity of air) 


Property 

Gr/Epoxy 

PZT-4 

Elastic Properties: 



Eu (GPa) 

132.38 

81.3 

E* (GPa) 

10.76 

81.3 

E,, (GPa) 

10.76 

64.5 

G* (GPa) 

3.61 

25.6 

G 13 (GPa) 

5.65 

25.6 

G n (GPa) 

5.65 

30.6 

V 12 

0.24 

0.33 

V U 

0.24 

0.43 

V23 

0.49 

0.43 

Piezoelectric coefficients (Iff 12 mV): 



d 3X 

0 

-122 

d 32 

0 

- 122 . 

das 

0 

285. 

^24 

0 

0 . 

Electric Permittivity: 




3.5 

1475. 


3.0 

1475. 

£ 33/^0 

3.0 

1300. 

Mass Density p(kg/m 3 ) 

1578. 

7600. 
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Table 2. Predicted Natural Frequencies of [p/0/90/0/p] square simply supported Plate. 

(C)- closed circuit; (O)- open circuit 


Mesh Density 

Discrete Layers 
(N-l) 

Fundamental Frequency, / ; c?/h p m 
10 3 Hz (kg/m) 1/2 



adi=4 (C) 

a/h=4 (O) 

Cdh=50 (C) 

a/h=50 (O) 

w-constant: 






4x4 

3 

148.154 

153.720 

249.860 

272.860 

8x8 

3 

145.323 

151.222 

236.833 

259.173 

12x12 

3 

144.768 

150.727 

234.533 

256.765 







8x8 

7 

142.469 

148.351 

236.785 

259.103 

8x8 

20 

142.221 

148.177 

236.784 

259.102 







w-variable: 






4x4 

3 

149.310 

154.591 

253.670 

276.185 

8x8 

3 

146.269 

151.964 

239.628 

261.703 

12x12 

3 

145.675 

151.531 

237.158 

259.655 







8x8 

7 

142.942 

148.692 

237.383 

259.630 

8x8 

20 

142.630 

148.586 

237.109 

259.895 







Exact: 


145.339 

145.377 

245.941 

245.942 
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(a) Smart Laminate ( b ) Discrete-layer laminate theory assumptions 


Figure 1 .— Typical piezoelectric laminate configuration . (a) Concept ; (b) Assumed through-the- 
thickness displacement and electric potential fields. 
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(c) Normalized Thickness (z/h) 


Figure 2. — F undamental mode in a thin [p/0/90/0/p] amply supported square plate (thickness aspect 
ratio a/h= 50). (a) Lateral displacement w/w^ (b) in-plane displacements u/w^ and 
v/w ^ ; (c) normalized electric potential . 
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Figure 3. — Fundamental mode in a thick [p/0/90/0/p] simply supported square plate (thickness aspect 
ratio a/h=4). (a) Lateral displacement w/w^ (b) in-plane displacements u/w^ and 
v/Wjuj ; (c) normalized electric potential ■ 
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Normalized Stress Normalized Stress 
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(b) Thickness (z/h) 

Figure 4. Modal stress fields of the fundamental mode in a thin [p/0/90/0/p] simply supported 

square plate (a/h= 50). (a) In-plane stresses o x / , <* y / (<* y )m«x » CT xy t ( a y)m*xl 0 5 ) 

interlaminar stresses / (<T y ) IM3t > / (a y ).». 
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Figure 5. — Modal stress fields of the fundamental mode in a thick [p/0/90/0/p] simply supported 
square plate (ot/h= 4). (a) In-plane stresses a x / (Oy)^ , c y / (a y / (Oy),,^ (b) 

interlaminar stresses / (g y ~) / (g y ),„. 
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Open Circuit 

•••• 

Closed Circuit 



Open Circuit 
Closed Circuit 


Figure 6. — Effect of electric boundary conditions on the fundamental mode of the thick [p/0/90/0/p] 
plate (<x/h=4). (a) in-plane displacements u/w m „ and v/w,^ ; (b) normalized electric 
potential ; (c) in-plane stress a x / (Oy)^ (d) interlaminar stresses /(Gy)^ G„ 

)max • 
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Figure 6. — (Continued) 
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(b) 

Figure 7. — Schematic configuration of cantilever plate beams with distributed piezoelectric patches. 

(a) [0/±45], plate with 15 piezoelectric actuators/side; (b) [0/90/±45], candidate plate with 
8 piezoelectric actuators/side. Dimensions in mm. 
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Figure 8. — Predicted induced deformations of the [0/±45], cantilever plate beam. Measurements by 
Crawley and Lazarus (1989). 
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Sensory Voltage (Volts) ^ Y ( m ) 
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(b) 

Figure 9. — First electromechanical mode of the [0/90/±45], cantilever plate with discrete PZT-4 
patches (1st bending), (a) Lateral deflections w/w,„; (b) modal sensory voltage (upper 
surface). 
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Figure 10 — Second electromechanical mode of [0/90/±45], cantilever plate with discrete PZT-4 

patches (1st torsion), (a) Lateral deflections w/w^; (b) modal sensory voltage (upper 

surface). 
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